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1) INTRODUCTION AND OVERVIEW 

A research proposal "STUDY OF EXPERIMENTS ON CONDENSATION OF 
NITROGEN BY HOMOGENEOUS NUCLEATION AT STATES MODELLING THOSE OF THE 
NATIONAL TRANSONIC FACILITY" was submitted to NASA-Langley dated 1 
January 1979 to cover a period of 1 April to 21 December 1979. Funding 
was received on 4 June 1979 to cover a period of 15 May 1979 to 14 
February 1980 with a later extension without additional funds to 30 
April 1980. The final report at hand covers the work done under the 
listed grant. For past history and motivation we refer to the text and 

references given in the Proposal. A cryogenic wind tunnel is based on 

•j* 

the twofold idea of lowering drive power and increasing Reynolds number 
by operating with nitrogen near its boiling point. There are two pos- 
sible types of condensation problems involved in this mode of wind tvmel 
operation. They concern the expansion from the nozzle supply to the test 
section at relatively low cooling rates, and secondly the expansion around 
models in the test section. This secondary expansion involves higher cool- 
ing rates and shorter time scales. In addition to these two condensation 
problems it is not certain what purity of nitrogen can be achieved in a 
large facility. Therefore one cannot rule out condensation processes other 
than those of homogeneous nucleation. In sum, the data on ^-condensation 
of previous work shown in Figure 1 (for references see the Proposal) in 


t E.C. Polhamus, R.A. Kilgore, J.B. Adcock and E.J. Ray, "The Langley 
Cryogenic High Reynolds Number Wind-Tunnel Program", Astron. & Aeron., 
October 1974. 
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conjunction with two typical isentropes may well be insufficient for 
accurate prediction of operating conditions for the Langley facility. 
To repeat in part from the Proposal. 

(a) Expected supercooling at higher pressures is not well 
known (supply pressure of 1 to 10 atm for a Mach number 
range of ~0.7 < M < 2 are anticipated for the transonic 
facility) . 

(b) The experimentally available range of the cooling rate 
effect on supercooling is not sufficiently wide to cover 
the needs of the transonic facility. Expansion rate 
effects in the free stream and the higher rates around 
models and their components must be understood. 

(c) The effects of purity on condensation of the nitrogen 
used for tunnel operation is not understood. Hetero- 
geneous and binary (or heteromolecular) effects must be 
studied to understand nucleation in a practical situa- 
tion for nitrogen available in large amounts . 

The question of what impurity level can be tolerated may well 
prove to be one overriding importance for the operation of the tran- 
sonic wind tunnel facility in the supersaturated regime. There is 
clear evidence that nitrogen expansions can be nucleated hetero- 
geneously resulting in little or insignificant supercooling. Hetero- 
geneous nucleation has been demonstrated in a qualitative manner for 

■f 

other gases and in a semi-quantitative way here at Yale for the 


t For references see the Proposal. 
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condensation of ethyl alcohol and also for sulfur hexafluoride, both 
present In non-condensing carrier gases. 

The study performed under this grant addressed Itself to some of 
the points raised. The summer of 1979 was used to compute isentropic 
expansions of nitrogen based on real gas (equilibrium) properties. The 
Berthelot equation of state was chosen as the thermal equation of state, 
and It was Integrated together with caloric data on N 2 at zero pressure 
taken from the literature. This procedure yielded the real gas Isentropes. 
The Berthelot equation Is found to be a substantial Improvement over van 
der Waal’s work in predicting the state of a real gas in the range of pres- 
sures of interest here. The more complicated equations such as that 
proposed by Beattie and Bridgman which are useful at much higher pressures 
can thus be avoided. The results of this calculation are given in the 
Appendix. It became clear that real gas effects are in fact negligible 
in the range of interest given, provided no condensation appears. 

In parallel, a study was carried out to determine design parameters 
of a small wind tunnel operated much like that used for previous conden- 
sation studies on sulfur hexafluoride in argon and argon in helium. 

A few of the results of more extensive studies are discussed in Section 
3 of this final report. An experimental setup such as shown schematically 


t B.J.C. Wu, P.P. Wegener and G.D. Stein, J. Chem. Phys. 68, 308 (1978). 
B.J.C. Wu, and G.A. Laguna, J. Chem. Phys. 71, 2991 (1979). 

ft B.J.C. Wu, P.P. Wegener and G.D. Stein, J. Chem. Phys. 6£, 1776 (1978). 
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in Figure 2 is envisioned. The gas supoly consists of bottled N 2 which 
is passed via a sonic orifice valve and additional control valves through 
a heat exchanger embedded in an insulated vessel filled with liquid nitro- 
gen. In this fashion supply temperatures down to about 100 K can be 
achieved. The once-through system operates continuously via a plenum 
chamber and two existing Stokes Microvac Model 412-H pumps operating in 
parallel at a combined displacement of 600 CFM. Instrumentation and 
methods of observing condensation are in principle along the lines of 
those to be discussed in the next Section of this report. 

The study revealed that it would be difficult to duplicate the low 
cooling rate of the Langley Transonic Facility as measured from the sup- 
ply to the test section. However, a wide range of cooling rates as 
typically to be expected around models can readily be produced in a 
range of nozzles operated to duplicate the thermodynamic states of the 
full-scale facility. The points raised initially in this Section could 
then be answered in a systematic study. 

Finally, a few experiments on condensation of nitrogen using existing 

tt 

equipment were suggested in the proposal for late summer of 1979. However, 
in our continuing contacts with Dr. Robert M. Hall of NASA-Langley — the 
technical monitor of our work — it became apparent that changes in the 
original plans might be useful as NASA's interests change. For one, the 

t In a separate study at our laboratory for the NASA Ames Laboratory it 
becomes apparent that our theoretical understanding of binary, hetero- 
molecular and heterogeneous nucleation is far from satisfactory. This 
fact suggests an empirical approach to solve practical engineering pro- 
blems at this time. 

tt These experiments cover the lowest cooling rates ever achieved in steady 
supersonic flow used for condensation studies of nitrogen. 
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primary condensation interest appeared increasingly to focus on the 
expansion around models. No chance could be taken to operate the test 
section at a supersaturated state. Therefore it was decided to operate 
at Yale as part of the subject grant a more extensive series of nitrogen 
condensation experiments. In turn consumables (primarily gaseous and 
liquid nitrogen), small parts, etc. were used in excess of the budgeted 
amounts. Moreover more manpower than originally planned was needed to 
prepare the condensation setup. The actual expenditures exceeded the 
budget and they will be submitted separately. 

The results of these experiments covered at least a part of the 
range of interest to the Langley facility. In particular the cooling 
rates were somewhat lower than those of previous experiments (Figure 1) 
and they resulted in a smaller supersaturation before condensation set 

4 » 

in. Moreover, it is believed that our method of detection of the onset 
of condensation by light scattering is more sensitive than the instru- 
mentations! approaches of previous work. Prior to our experiments nitro- 
gen condensation was observed by pressure measurements. Thus our data may 
lead to an earlier detection of phase change on a given isentrope. This 
work will be discussed in the next section. 

2) NITROGEN CONDENSATION EXPERIMENTS 
a) Experimental Results 

The experiments to be described next covered the range of states in 
the vicinity of Isentrope # 2 in Figure 1. The cooling rates were of the 
^.rder of 0.1 C/us. Apparatus close to that shown in Figure 2 was used, 
however the light beam of a Helium/Cadmium laser (442 nm) was passed 

t See footnote on previous page. This statement applies to the pressure 
range of our experiments. 


beyond the nozzle exit via a first surface mirror and a prism at right 
angles to the flow direction. In turn the photomultiplier was mounted 
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at 90 to the Incident light beam perpendicular to the flow direction. 

Bottled nitrogen of purity 99.9982 was passed via a sonic orifice and a 

control valve through three separate copper coils embedded in an Insulated 

heat exchanger filled with liquid nitrogen. ^ The cooled nitrogen gas 

entered the supply section as shown in Figure 2 and the double walled 

supply was cooled with liquid nitrogen. The arrangement is shown in 

Figure 3 with details of the laser seen in Figure 4. The supply tempera- 

0 

ture was measured by a thermocouple versus a cold junction at 0 C. A 
Wallace and Tieman aneroid pressure gauge in parallel with a Baratron 
were used to measure the supply pressure. Previously it had been shown 
that the supply temperature remains constant across the supply chamber 
of the wind tunnel. The range of temperatures covered was 107 < T q < 124° K. 
Higher temperatures were difficult to achieve because individual coils 
could not be removed from the circuit without extensive redesign. 

The nozzle used was designated as Nozzle B as listed in Table 1. It 
is moreover shown full scale in Figure 5. A free jet issued into an open 
chamber at the nozzle exit (Figure 2) and the jet could be observed through 
two windows. Several pressure taps connected to a single tube are located 
3 mm \ipstream from the nozzle exit (Figure 5). This pressure is designated 
as p e and it is measured with a Baratron pressure gauge which was calibra- 
ted with a Wallace and Tieman mercury precision manometer. Finally, the 
pressure in the plenum chamber was kept at or below the nozzle exit pres- 
sure. 

t The N 2 was Grade P.P. ("prepurif ied") with max. moisture 3 parts/million, 
max. O 2 5 parts/million (Linde Corp.). 




The nozzle was first calibrated by measuring the exit pressure as 
a function of the supply pressure for a room temperature supply. The 
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results of this calibration are shown in Figure 6 together with the com- 
puted effective area ratio. We note the Reynolds number effect with 
decreasing exit pressure ratio (i.e. higher Mach numbers) at higher sup- 
ply pressures. At p Q *■ 1 atm, the Reynolds number per unit length at 
the nozzle exit is about Re/£ * 10 7 m” 1 at M * 2.5. Unfortunately a 
similar calibration does not exist for nitrogen with a "cold" supply 
because of a malfunction of the pressure connection that was subjected to 
the low temperatures. 

The condensation experiments were conducted with continuous flow. 

Test section and heat exchanger were precooled by passing ^ low flow 
rate through the apparatus. The supply temperature was varied by alter- 
ing the mass flow rate in the heat exchanger by a variation of the supply 
pressure.^ After the desired temperature was obtained, the supply pres- 
sure was raised until a light scattering response was barely perceived 
above the signal given by the background light and the electronic noise. 
Lowering the pressure led again to a disappearance of light scattering 
and repeated pressure changes up and down provided a check of the opera- 
tion. The voltage from the photomultiplier, the supply pressure, Baratron 
output and the supply temperature were recorded on Hewlett-Packard strip 
chart recorders and simultaneously read on digital voltmeters. 

t Independent temperature control by separate operation of parts of the 
heat exchanger will be essential in future experiments. 
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A series of condensation experiments as described were performed 
from 28 November to 2 December 1979. Of this series 17 reliable data 
points on the onset of condensation of nitrogen in Nozzle B (Table 1, 

Figures 5 and 6) were obtained. The raw daca and derived results of 
these experiments are listed in Table 2 as explained in the legend of 
Table 3. The data reduction method will also become apparent from 
Table 3. The main points of the procedure are as follows. For each 
experiment column (1) at a measured value of steady flow with p Q (2) , 
and T o (3), (4), the onset of condensation was noted at the nozzle exit 
by a barely perceptible light scattering increase above the background 
signal. Because of the unfortunate absence of a "cold" calibration, 
two limiting values of the corresponding onset state of condensation 
were determined from the raw data. One limit is given by assuming the 
absence of viscous effects, i.e. flow without a boundary layer. The 
other limit is obtained by relying on the room temperature nozzle cali- 
bration of Figure 6. Clearly the actual value on a given Isentrope 
will be found somewhere in between the two extremes (note later dis- 
cussion). The onset of condensation in these experiments at M * 2.3 does nut 
lead to thermal choking permitting the straightforward calculation of 
states on an isentrope without regard to heat addition by condensation 
whose value is much below the critical oue. Finally, all results of 
Table 2 were computed by treating nitrogen as a thermally and calorically 
perfect gas.* 

The experimental results of Table 2 are given in Figure 7 . The 

+ 

The justification for this piocedure is discussed in the Appendix. 
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saturation vapor pressure, p^CT) with respect to a flat surface is also 
plotted. Seven sample isentropes are shown for supply pressures of 0.5 
to 9 atm at different supply temperatures. Mach numbers are marked off 
on these isentropes. Experiment #32 at - 129K and p Q * 1.5 atm is 
also indicated by an isentrope with the limiting condensation states 
markea oft. All other experiments are indicated by their number at the 
supply state and a short isentrope that brackets the onset of condensa- 
tion within the two assumptions previously given. The mean value of 
the limiting experimental results is shown by a dashed line which indi- 
cates an adiabatic supercooling^ of roughly 10 *C shifted in parallel to 
the vapor pressure curve. 

The cooling rate will next be calculated. With relatively small 
changes of the Mach number, flow speed, etc. between saturation and con- 
densation, it is reasonable to define an average speed between these 
two states. The state at saturation can be roughly located in the noz- 
zle based on thi: effective area distribution A/A* » f(x) from the previous 
nozzle calibrations. The onset of condensation is, of course, found 
at the nozzle exit in view of our experimental method. Choosing the mean 
value curve of Figure 7 for the onset of condensation, the results shown 
in the upper half of Table 4 emerge. We note that the adiabat super- 
cooling, AT^, has the narrow range of 

t National Bureau of Standards, Circular 564, (1955). 

tt The adiabatic supercooling, AT ., is defined as the temperature dif- 
ference between saturation and condensation measured on an isentrope. 

ftf These calibrations were recomputed since they had been carried out 
for y ■ 1*67 rather than y m 1.40 (see references in the footnote 
on p.3). The effective area ratio at the exit determined from the 
present calibration with N 2 agreed with the previous calibration within 
5% at p • 100 torr and 1.5% at p - 800 torr. 
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10 to 11 *C for the 4 exper intents evaluated. Clearly these results 
are no better than the initial estimates of the onset state and we con- 
clude AT^ * 10 *C. 

b) Discussion of Results 

The following discus*-' ->n has two aims. For one, we shall view our 
nitrogen experiment in the context of previous work. On the other hand, 
the experiments possibly to be pursued in the future will be subject to 
the same methods to be discussed here. Essentially similar work in a 
much extended range and with high accuracy can be expected. At first 
some remarks on the determination of the onset of condensation by light 
scattering. 1 No quantitative calibration was carried out. Owing 

to careful masking, use of a dark room and application of optical 
black on many components, the background light was minimal. Incipient 
condensation led to a large increase in the Rayleigh light scattering 
signal. The photomultiplier tube used, the tube voltage, the light 
scattering geometry, the He/Cd laser, and the thermodynamic state at the 
nozzle exit were all in a range similar to that of several previous 
investigations. Consequently we expect the measured increase of the scat- 
tered light to reflect clustering with concentrations of the order of 
10 11 clusters/cm 3 with corresponding droplet radii of the order of 15 to 
20 Angstroms. The mass fraction condensed at that location is not 

t This discussion leaves out the unknown exact state as caused by the 
absence of a "cold" nozzle calibration. (See Tables 2 and 3 and 
Figure 7). 
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readily determined, in fact this is the type of result which can later 
be estimated more accurately. However, with the above assumption and a 
typical density value (e.g. Exp. it 32, Yale, Table 4) at the nozzle exit, 
the mass fraction of condensate, g, that can be detected is estimated to 
be 10“ 6 < g < 10” 5 . This value is at least two orders of magnitude 

lower than could be observed by a deviation of pressure from an isen- 
tropic expansion. As stated before, we may therefore find more realistic 
values for the adiabatic supercooling than noted in the past with a 
corresponding lower value seen by us. 

It is thought, however, that the adiabatic supercooling may be 
slightly larger than shown in Table 4. This is a consequence of the 
fact that the boundary displacement thickness of the cooled nozzle may 
be somewhat lower than the one predicted by the mean values of Table 4 
as shown in Figure 7. It is further interesting to note that the super- 
cooling of Table 4 roughly agrees with that found at our laboratory in 
1978 for pure argon using a similar apparatus. The surface tension of 
liquid argon (T * - 188C, a * li erg /cm 2 ) is in fact not too different 
from that of liquid nitrogen (T * - 203C, c * 11 erg/cm 2 ). The vapor 
pressures and indices of refraction are not too different. Consequently 
onset conditions of condensation may not be too different yielding 
similar differential cross sections. 

These thoughts finally bring us to a comparison of our results — 
always based on the mean value as indicated by the dashed line in Figure 
7 or the sample values of Table 4 — with previous work. In Figure 1 
taken from the Proposal we show nozzle and free jet data. The nozzle 
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data refer to the work at the Cal. Inst, of Tech, in 1952 (Arthur, 

Buhler, Willmarth and Nagamatsu, see Proposal) and the work of the 
Applied Physics Lab. of JHU in 1952 (Faro et al., see Proposal). The 
free jet data came from recent (1978) German experiments (Dankert and 
Koppenwallner , see Proposal). The curve marked "Goglia" in Figure 7 
represents his Interpolation of his data and those of the CalTech group. 
However, the Goglia data^ at low Mach numbers but higher pressures are 
most suitable for comparison with our results. We note that for conden- 
sation onset pressures below about 200 torr — roughly where the word 
mean appears in Figure 7 — Goglia connects his data points with those 
of the CalTech group. Shifting our mean curve slightly to the left, 
i.e., uhe cold nozzle site we find that an extension to higher pressures 
shows reasonable agreement of our data with Goglia" s values. In fact an 
estimate of the adiabatic supercooling based on three of Goglia" s experi- 
ments at higher pressures and temperatures given in Table 4 shows roughly 
the same adiabatic supercooling that we found. This comparison brings 

us to the final discussion of the cooling rate the second empirical para- 

tt 

meter next to pressure level, that governs the adiabatic supercooling. 

In general we have drastically different cooling rates for different 
experimental methods of study of nucleation. The expected ranges are 
roughly as follows: 

t In addition to the readily accessible reference given in point (1) of 
Table 5 note: Gennaro Louis Goglia, "Limit of supersaturation of nitro- 
gen vapor expanding in a nozzle," Doctoral Dissertation, University of 
Michigan 1959. 

ft We use the adiabatic supercooling since small errors in pressure mea- 
surement cause large errors in the calculation of the supersaturation. 
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Estimated 

Metho d Cooling Rate Range (C/ps) 


atmospheric motion 

10" ; 10 

to 

10’ 7 

diffusion cloud chamber 

10" 8 

to 

10’ 7 

piston cloud chamber 

10’ 7 

to 

10’ 5 

shock tube (expansion fan) 

1 

O 

H 

to 

10’ ' 1 

wind tunnel^ 

10’ 2 

to 

10 1 

shock tube (tailored interface) 

- 

0 

10 


molecular beam 

10 2 

to 

10 4 


In this list we note the extreme rates of cooling in free jets of 
molecular beams, consequently large supercooling is expected as observed 
in Figure 1. The shaded range in Figure 1 and 7 includes the results of 
the diverse sources mentioned. Most of the data were taken in hypersonic 
tunnels hence the overlap of Goglia's extrapolation at lower pressures. 
Our wind tunnel experiments were in contrast performed at cooling rates 
of about 0.10 C/ms as computed in Table 4 for four experiments. We see 
in the table that saturation takes place at M - 2, hence with condensa- 
tion at M - 2.6, a lower cooling rate arises. It is this feature of the 
precooled, low Mach number experimental apparatus that can be varied to 
produce a wide range of cooling rates applicable to the flow fields 
about models in the transonic facility. Goglia's cooling rates, also 
estimated from his data in Table 4 are abouc twice ours. In line with 


t Supply to test section, includes hypersonic tunnels. 
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past experience, a factor of two of this variable is not too important; 
thus no significant effect on the adiabatic supercooling arises between 
the two series of experiments. However, more accurate and extended 
experiments are clearly needed. In sum, our condensation data 

provide results of interest in the lower pressure range of the 
transonic facility, 

have detected the onset of condensation at much smaller mass 
fractions than noted previously, and 
demonstrated that an extension of the techniques used could 
cover the entire range noted in Figure 1 including a 
substantial variation cf the cooling rate. 

3) STUDY OF PROPOSED CONDENSATION EXPERIMENTS. 

A study is proposed to be carried out at Yale University where 
nitrogen expansions of the NASA Langley National Transonic Facility shall 
be duplicated in a small continuously operating wind tunnel covering a 
small fraction of the cross-sectional area of the prototype. Experiments 
are to be performed covering the range of supply states, Mach numbers and 
cooling rates around models of the full scale transonic tunnel to observe 
if supersaturation of the nitrogen can conservatively be expected. Pre- 
liminary experiments of this kind carried out with similar equipment were 
described in the previous Section 2 of this Final Report. 

It is proposed to vary the conditions important to condensation in 
these experiments, by again (TABLE 6) varying the test gases at first from 
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ultrapurif ied nitrogen to prepurified nitrogen. 1 In addition it will be 
important to add systematic studies involving binary and heterogeneous 
nucleation. In such experiments small metered amounts of H 2 O and/or aero- 
sols of known geometry will be added to the nitrogen. In all cases the 
supply section, however, will have pure vapors (no dropletsl) present. T 

The proposed test facility will be instrumented by a combination of 
pressure measurements, Rayleigh light scattering experiments (light beams 
across the flow and on the centerline) , schlieren photography and possibly 
interferometry.^ The experimental results will be compared with predic- 
tions of condensation based on solutions of the equations of motion in 
conjunction with the theory of nucleation and appropriate droplet growth 
laws. 


The System 

A wind tunnel of the kind shown in Figure 2 would be designed and 
constructed. Its basic features are discussed in Section 2 and the 
references given therein. The maximum size of the tunnel is limited by 

the capacity of our two existing Stokes Microvac Model 412-H pumps opera- 

ffi' 

ting in parallel. Their joint capacity for 0.01 < P inta j ce < * atin is 
2 x 270 CFM - 540 CFM - 0.25 m 3 /s. Assuming a range of initial conditions 
in the nozzle supply as shown by two limiting isentropes in Figure 1, the 


t No significant differences are to be expected between both purified 
test gases as noted before in previous investigations of similar pro- 
blems. See footnotes on p.3, also Wegener et al., Phys. Fluids 15, 

1869 (1972) and Wegener and Wu, Faraday Discussion of the Chemical Soc., 
No. 61 (1976) in which binary nucleation experiments in supersonic flow 
are discussed. 

ft We will consider the use of our Zeiss Mach-Zehnder Interferometer. 

ttt The capacity of our vacuum system can, of course, be increased by 
additional pumps if desirable. 
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data of TABLE 6 were computed. Here It is conservatively assumed that 
the throat area equals 1 cm 2 and that the pump intake pressure equals 
the test section pressure while the gas is isobarically heated to room 
temperature in the return pipes leading to the vacuum pumps below the 
floor of the laboratory. In other words we assume no pressure recovery 
in a diffuser. Note that the proposed tunnel is appreciably larger 
(compare Figure 5) than that used in the experiments described in Section 
3. The estimated N 2 gas consumption and its cost is also shown in TABLE 

6. "*"*" A single experiment, i.e. the accurate determination of one point 
of incipient condensation in a p(T) -diagram such as that shown in Figure 

7, takes no more than about five minutes running time. 

The N 2 gas is precooled in a liquid nitrogen heat exchanger similar 
to the one currently in use (see bottom of Figure 3). However, the gas 
flow will be divided into a number of different copper coils (3 or 4) 
whose discharge rate may be controlled independently. Moreover an addi- 
tional gas feedline by-passes the heat exchanger, and all flows are mixed 
before entry into the tunnel supply. This mixing cakes place in an insu- 
lated length of pipe between heat exchanger and tunnel supply computed for 
full mixing by turbulent pipe flow. In this fashion it will be possible 
to control the supply temperature of the wind tunnel in a wide range as 
needed to cover the demanded isentropes. Heat transfer calculations based 
on existing property data and appropriate values of heat transfer coeffi- 
cients, thermal conductivity, etc. reveal a required total length, L, of 

t The cost figures are subject to price escalation. 

it *5 

ft It is expected that reliable experiments can be performed at A < L cm-. 
Recall that A and consequently cost are linear functions of the throat 
area. It may well be desirable to build a second nozzle with a smaller 
throat for the high pressure isentrope studies. 
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copper tubing (ID - 1 cm), fully immersed in an Insulated box filled with 

O 

liquid nitrogen at atmospheric pressure (~ 77 K),of approximately 5 < L < 8m. 
The total length depends on supply conditions and the choice of separately 
operable coils and a by-pass will permit adjustment to a given supply tem- 
perature. These values are in line with our past experience (See Section 3). + 

A number of substantial differences in design of the test section in 
comparison with the existing system (Figures 2 to 4) will be made. For one, 
the supply section will contain a fixed thermocouple rake to avoid survey- 
ing the supply temperature uniformity with a moveable thermocouple. Windows 
to observe the prescribed absence of condensation determined from light scat- 
tering in the supply will be available. All windows will be connected to 
the supply and test sections by evacuated chambers in order to avoid fogging. 
The entire supply section will again be double-walled to be cooled with 
liquid N 2 before and during actual tunnel operation. The nozzle will be of 
the type shown in Figure 5 with a continuous expansion. A cubic throat 
curve will be joined to a plane diverging wall to achieve a continuous 
second derivative of the surface to insure a smooth turbulent boundary 
layer growth. The nozzle will be equipped with static pressure taps to 
calibrate the flow without resorting to a nozzle calibration in a separate 
setup outside the cooled test section as done in the past. This calibra- 
tion will be monitored at all operating conditions. The nozzle exit will 
extend into a large chamber. That chamber will have circular windows on 
all sides far removed from the jet to observe the jet visually and by the 

t In fact the practical difficulty in the operation of such heat exchangers 
(or precoolers) lies in achieving temperatures appreciably above those of 
liquid nitrogen at its boiling point. 
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optical techniques mentioned before. In this fashion fogging and other 
low temperature problems can be minimized. Moreover a number of fiber 
optic channels will be incorporated in the test section for attachment 
of photomultipliers on the outside at room temperature. In this fashion 
the onset of light scattering in the nozzle itself can be monitored. 
Currently this point was shifted to the nozzle exit by manipulation of 
the supply state. A more powerful blue laser than used in the previous 
work would also be useful (e.g. Liconix Helium Cadmium Laser Model 4050 
giving 40 mW minimum at 442 run) . 

The Flow 

For the system described and the supply conditions of TABLE 6, the 
Reynolds numbers per unit length at M ■ 2 in the free stream will have 
the approximate range of 10 4 <Re/£ < 5*10 5 nr 1 . Estimates of the dis- 
placement thickness of turbulent boundary layers at these conditions at 
a typical location about 3.0 cm from the throat lead to effective cross- 
sectional areas roughly 95% to 98% of the geometrical area for M - 1.4. 
This is a more favorable result than found from the calibration with a 
"warm" tunnel shown in Figure 6. At any rate, for each condensation 
experiment a calibration of the nozzle at identical supply conditions 
will be performed. 

Finally, the cooling rate must be discussed. For simplicity, and 
in contrast to the results of TABLE 4 with saturation achieved at M > 1 
Csee Figure 7) , we choose the throat region at that part of the nozzle 
giving the highest cooling rate. 


t The laboratory permits operation in total darkness. 
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At the nozzle throat denoted by an asterisk we have 


d(T/T Q ) ]* 
dt 


(R*h*)“* 


( ^ 1)a o 

[<y+l>/2] 2 


» 


with 

R* radius of curvature at the nozzle throat, 
h height of throat, 

supply sound speed, and 
y ratio of the specific heats. 


For nitrogen and with y « 1.A0, the above equation may be simplified to 
give 

(dT/dt)* - - (R*h*)”* T 3/2 5.67 (°C/s) , 

o 

with T q in K and the length R* and h* in meters. 

For convenience we can approximate the throat and the following 
geometrical expansion in the nozzle by an arc. Using the equations given 
we can use typical throat radii (e.g. R* * 1 m) , nozzle lengths and nozzle 
heights compatible with the system described before. Cooling rates near the 
throat of 0.01<-dT/dt<0.1°C/s can be achieved. These values are even lower 
than those 'hown in TABLE 4 for the preliminary experiments. As one 
example we choose a nozzle with R* * 1 m and h* * 10” 2 m (for a 1 cm 2 throat 

O O 

area), and we assume T q - 100 K. Our equation gives -dT/dt = 0.06 C/s. 


t See p. 707, Peter P. Wegener "CONDENSATION PHENOMENA IN NOZZLES", 
Progress in Astronautics and Aeronautics, pp. 701-724, Academic 
Press, Inc., New York (1964). 
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Nozzle 



TABLE 1 ; Geometry of various nozzles used in condensation 
experiments at Yale. Nozzle B used for pure 
nitrogen condensation experiments reported here. 
Similar nozzles would be proposed for extensive 
studies. 
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TABLE 3: 


Legend fcr TABLE 2. + 


(1) Number of experiment. 

(2) Supply pressure , p Q . 

(3) , (4) Sutoly temperature, T . 

0 

(5) p £ refers to pressure at the exit of the nozzle (Jt«8 cm, TAiiLE 1). 

P e /p Q as shown was measured for pure nitrogen (no condensetlon) as 

a function of p Q at a supply temperature of roughly 290 K ( see 

Figure 1) . These are che values given in column (5) . Boundary layer 

effects for this "hot" case — (compare T -values of columns (3) and 

0 

(4) — lead to effective nozzle exit area ratios (see Figure 1} that 
are smaller than those of the low temperature case. For low T - 
velues no calibration exists. 

(6) p found by multiplying values of columns (2) and (3) . 

* (y-1) /y 

(7) T , estimated exit tempereture, found from T - T (p /p ) ' 

€ ^ O ® O 

with T from column (3) and p /p from column (5) with y ■ 1.4. 
o e o 

The values in columns (6) and (7) give one limiting condition of the 
thermodynamic state of expected condensation onset on the isentropc. 

(8) p g computed for inviscid flow. From TABLE 1 we find the exit Mach 
number of Nozzle B as M • 2.75, Consequently p./p “ 0.0398 for 

y ■ 1.4. Multiplying this fixed value with those of column (2) gives 
the values of column (8). 

(9) T computed like the T of column (7), however with fixed p /p * 0.0398. 

c €5 f» O 

Columns (8) and (9) consequently give a second limit for the condensa- 
tion onset conditions on an isentrope, the zero boundary layer limit. 

(10) M derived from column (5) . The corresponding Mach number for Inviscid 

flow as derived from the area ratio is given by - 2.75. 


t Condensation results of TABLE 2 are shown in Figure 7. 

ft It was separately shown that the flow car. be treated as that of a calori- 
cally and thermally perfect gas. 


Saturation Onset (mean) 



T ABLE A: Thermodynamic state at saturation and Bean values 

at the onset of condensation for sample experiments 
of this work and that of G.L. Goglia. 


TABLE 5: 


Legend for TABLE 4. 

(1) Number of experiment. For supply state, etc. see TABLE 2 for the 
Yale experiments. For the Goglia experiments see e.g. G.L. Goglia 

and G.J. Van Wylen, Trans. ASME 83, Ser.C, 27, (J. Keat Transfer), 1961. 

(2) , (3), (4) Thermodynamic state at saturation based on p oo (T) of NBS, 

Circular 564, (1955). 

(5), (6), (7) Thermodynamic str.e at the onset of condensation. Mean 
value of states for inviscid flow and "hot" calibration (see TABLE 
3). 

(3) Average flow speed between saturation and the onset of condensation 
u(m/s) . 

(9) Adiabatic supercooling AT a ^(C), i.e. temperature difference between 
saturation and condensation onset. 

(10) Average cooling rate dT/dt (C/ys) computed by dividing AT by the 
time of motion of a fluid element at speed u (8) between the loca- 
tion of saturation in the nozzle and the nozzle exit. Note that 
condensation is forced to appear at the nozzle exit. 
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TABLE 6 : Consumption rates of a wind tunnel operated with N 2 with a 

throat area of 1 cm 2 in # the range of isentropes of Figure 1. 
Discharge rates of N 2 (m) at the supply conditions indicated 
and use of liquid N 2 to precool the gas to the supply state 
are shown. Note additional cost estimates. 


NITROGEN CONDENSATION 



FIGURE 1: Results on pure N 2 condensation in nozzles 

and free jets from several sources (see 
proposal) . Isentropes 1 and 2 bracket the 
conditions for operation of the National 
Transonic Facility. 






« 


setup used to study condensation. 

A similar arrangement would be proposed 
for an intensive study. 



FIGURE 3 ; Wind tunnel (flow from 

left to right) and plenum 
chamber. Note liquid N 2 
cooler at bottom. 


FIGURE 4 


He-Cd laser (blue) and 
optical arrangement to 
determine light scat- 
tering at nozzle exit. 
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FIGURE 5 : Nozzle B (also note Table 1); Full scale. 

A/ A* * 3.3 at nozzle exit. All dimensions 
in centimeters. 
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APPENDIX: Isentroplc Expansion of Nitrogen using Berthelot's 
Equation of State . 

List of Symbols 

A area 

c constant 

0 ^ specific heat at constant pressure 
c ideal gas (thermally ideal) specific heat, Eq. (9) 

h enthalpy 

h Q enthalpy at reservoir conditions 
p pressure 

P c critical pressure 
p density 

R gas constant 

s entropy 

T temperature 

T c critical temperature 
u flow speed 

v specific volume 

x distance 
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Summary 

The steady, one-dimensional, isantropic expansion 
of nitrogen is formulated using the equation of state 
suggested by Berthelot. The gasdynamic variables for 
supersonic nozzle flow are calculated assuming that no 
condensation takes place even though the nitrogen may 
be in supersaturated states. The results show the 
Berthelot gas model to be close to (less than 0,2% 
deviation) both the real gas calculations of Adcock 
(Ref. 1) and the ideal gas calculations. Therefore, 
the ideal gas equation may be used with sufficient 
accuracy in the temperature and pressure range of 
interest. 

Derivation 

The conservation equations of steady state, one-dimensional, 
cid nozzle flow take the form: 

continuity: pAu = m = constant. 

du . 1 dp A 
momentum: u + T *, ■ U. 

dx p dx 

1 o 

energy: + h = h ■ constant. 

4. O 

Combining Eq. (2) and the differential form of Eq. (3) gives 

— dp = dh, 

P 


:’nvis- 

( 1 ) 

( 2 ) 

(3) 

(A) 
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which implies ds * o by virtue of the thermodynamic relation 

dh - Tds H* 7 dp. 


(5) 


Hence, the momentum equation may be replaced by the isentropic con- 
ditions : 

s = s = constant. ( 6 ) 

o 


Thus, the solution of tne ret of equations (1), (3), and ( 6 ) requires 

expressions for the enthalpy, h(p,T), the entropy, s(p,T), and the 

thermal equation of state, o(p,T). In fact, the ideal gas specific 

heat c (T), and p(p,T) are all that is needed to generate the func- 
po 

tions h(p,T) and s(p,t). 

The Berthelot Thermal Equation of State*can be written 


pv = RT 


1 




(7) 


where the subscript c denotes the critical properties. For nitrogen 
we have 

T * 126.2 °K, and p ■ 33.5 atm. 
c c 


Thus, Eq. (7) can be written as 



* Ref. 3, p. 296. 


( 8 ) 
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where 


9 

128 


- 3.485x10 


-4 


K 

non 'Jg 


c, * 6 I 2 ■ 9. 556 x 10 4 (°K). 
i c 


The ideal gas specific heat capacity at constant pressure, 
from Jacobsen (Ref. 4) and it is expressed by 


c is taken 
po 


c 

po 


(T) 




(9) 


The 7 constants in Eq. (9) are tabulated by Jacobsen based on experi- 
mental results. The results obtained with Eq . (9) are plotted in Fig. A1 
in the form of a deviation of real from ideal heat capacity values given 
in percent. We note that for 30 <T- UO^K, nitrogen is essentially a 
calorically ideal gas. The ideal value of c /R " 3.5 is give: by in 
Eq. (9). 

Next the enthalpy and entropy are derived by following the procedure 
described in Ref. (2) on pp. 10-32 to 10-39. The path of the integration 
shows 3 steps. The results are 


h(p,t) 




( 10 ) 


and 


3 (p,T) * - Rlnp + R 2 -j TIT.Jj* + N 4 1"(T) 

L il«4 


- 2Rc^C2 



(11 


The speed of sound is required to determine the Mach number. In 
general for an equation of state p ■ p(p,T), the speed of sound, c, is 
given by 


c 



(12 


where 



(to\ (1 

1 

1 

* 

«l 

W/ T \ n ) v v 

py T v j 


(13 


For the Berthelot gas, the required derivatives in Eq. (13) are: 



c 2 rt 

2 ’ 
P 


and 






6c,c, 



(14 


(15 


(16 


(17 


- R 




Finally the sec of Eqs. (1), (3) and (6) can be solved using 
temperature as the independent variable. From a given initial state, 
(P o ,T q ), the et.tropy of the system, s q , is calculated using Eq. f 11) - 
To find the pressure at any other temperature, T, Eq. (11) is solved 
by trial and error. Once the pressure has been found, the remaining 
variables are found in a straightforward way. The computer was used 
to facilitate these calculations. 

The results are shown in Table Al. Less than 1% deviation in 
pressure and temperature from the ideal gas values are found over the 
range of interest, however, Adcock's real gas calculations (Ref. 1) 
also show less than 1% deviation from the ideal gas values. Hence, in 
the range of interest, the ideal gas equation may be used with suffi- 
cient accuracy. This result preatly simplifies the determination of 
the states of observed onset of condensation permitting the use of 


standard flow tables. 
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FIGURE Al : Deviation in 1 of the specific heat of nitrogen 

computed from Eq.(9) in relation to the ideal gas 
value. 
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0.9967 

0.3703 
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0.5519 


T/TO 


1 .OOCC 

C .909* 
C.95fet 

C .996* 

VM\l 

C .96 1* 
C .9 756 
C .960. 
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0 .9c* 12 

l .,*46 

3.6032 

157.18 
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0. cH-7 

1-5 . 7 9 

0.8*32 
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0.8 ?21 
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0.4309 
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0. ft* 31 
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0.9649 1.5873 
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0.92. 6 1.13*9 

0,4137 1.0439 

0 .®0e3 1.0621 

0.890* 1.0381 


MG >15.7010 CAL/Gh 

C P/PO T/TO y/VG A/A® 

AEIATIVE TO I0EAL GAS VALUES 
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0.7283 
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0.700* 
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0.5942 
0. 5590 
0. 59a9 
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0.99J.) 
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0 .5617 
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I.tOOO 0.9997 
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0.5Q95 C.9»98 

O .5995 C .9998 
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v . 56*8 
-. * -9 5 
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C .9998 
C. 3998 
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0.9969 
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